79-FM-26 
Rev. 1 


JSC-16029 


OPS MCC 
Level B/C 

Formulation Requirements 


Area Targets and Space Volumes Processor 


( N AS A-TM— 80839) OPS HCC LEVEL B/C N80-13840 

FORMULATION REQUIREMENTS: AREA TARGETS AND 

SPACE VOLUMES PROCESSOR (NASA) 140 p 

HC A 07/MF AO 1 CSCL 12A Onclas 

G3/64 42722 


Mission Planning and Analysis Division 
October 1979 


NASA 

National Aeronautics and 
Space Administration 

Lyndon B. Johnson Space Center 

Houston. Texas 


/c _ ' 




79-FM-26 


JSC- 16029 


SHUTTLE PROGRAM 

OPS MCC LEVEL B/C FORMULATION REQUIREMENTS 
AREA TARGETS AND SPACE VOLUMES PROCESSOR 


By Mack J. Bishop, Jr. 
Computer Sciences Corporation 


JSC Task Monitor: Alfred A. Menchaca f- A/ 

Flight Planning Branch 


Approved: 


Edgar C# Lineberry, Chief 
Mission Planning and Analysis Division 


Approved: 



Ronald L. Berry, Chief 
Mission Planning and Analysis Div n 



Mission Planning and Analysis Division 
National Aeronautics and Space Administration 
Lyndon B. Johnson Space Center 
Houston, Texas 
October 1979 



ACKNOWLEDGMENTS 


The author wishes to gratefully acknowledge the technical 
inspiration and guidance provided by A. A. Menchaca of JSC 
and Messieurs R. Rich and H. Norman of IBM. Their patience 
and ideas contributed greatly to the contents of this document. 


iii 


PRECEDING PAGE Bl,AN7{ NOT FILMED 



CONTENTS 


Section Page 

1.0 INTRODUCTION 1-1 

2.0 AREA TARGETS AND SPACE VOLUMES CHARACTERISTICS 

AND CONTAINMENT CRITERIA 2-1 

2.1 EARTH- REFERENCED CIRCLES 2-4 

2.1.1 Procedure 2-4 

2.1.2 Equations 2-4 

2.1.3 Assumptions and Limitations 2-11 

2.2 CELESTIAL CIRCLES 2-12 

2.2.1 Procedure 2-12 

2.2.2 Equations 2-12 

2.2.3 Assumptions anu Limitations 2-14 

2.3 EARTH- REFERENCED POLYGONS 2-15 

2.3.1 Procedure 2-15 

2.3.2 Equations 2-15 

2.3.3 Assumptions and Limitations 2-20 

2.4 CELESTIAL POLYGONS 2-21 

2.4.1 Procedure 2-21 

2.4.2 Equations 2-21 

2.4.3 Assumptions and Limitations 2-23 

2.5 EARTH- REFERENCED SPACE VOLUMES 2-25 

2.5.1 Procedure 2-25 

2.5.2 Equations 2-25 

2.5.3 Assumptions and Limitations 2-31 


v 


PRECEDING PAGE Bl.ANK NOT FILMED 



Section 


Page 


2.6 CELESTIAL-FIXED SPACE VOLUMES 2-32 

2.6.1 Procedure 2-32 

2.6.2 Equations 2-32 

2.6.3 Assumptions and Limitations 2-34 

2.7 SUMMARY OF EQUATIONS 2-35 

2.7.1 Earth-Referenced Circles 2-35 

2.7.2 Celestial Circles 2-36 

2.7.3 Earth-Referenced Polygons 2-36 

2.7.4 Celestial Polygons 2-38 

2.7.5 Earth-Referenced Space Volumes .... 2-39 

2.7.6 Celestial-Fixed Space Volumes .... 2-40 

3. 0 SEMI ANALYTICAL ALGORITHM TO PREDICT AOS AND 

LOS TIMES 3-1 

3.1 COMPUTING OUTER BOUNDARY CONE PARAMETERS. . . 3-4 

3.1.1 Earth-Referenced Circles 3-4 

3.1.2 Celestial Circles 3-5 

3.1.3 Earth-Referenced Polygons 3-5 

3.1.4 Celestial Polygons 3-7 

3.1.5 Earth-Referenced Space Volumes .... 3-7 

3.1.6 Celestial-Fixed Space Volumes 3-9 

3.2 DETERMINING WHETHER AOS AND LOS TIMES ARE 

POSSIBLE 3-11 

3.2.1 Celestial-Fixed Targets 3-11 

3.2.2 Earth- Referenced Targets 3-14 

3.3 PREDICTING AOS AND LOS TIMES 3-17 

3.3.1 Celestial-Fixed Targets 3-17 


vi 



Section Page 

3.3.2 Earth-Referenced Targets 3-21 

3.4 ASSUMPTIONS AND LIMITATIONS 3-29 

4.0 FUNCTIONAL OVERVIEW OF AOS AND LOS TIME 

COMPUTATIONS 4-1 

4.1 EARTH-REFERENCED AREA TARGETS AND SPACE 

VOLUMES 4-2 

4.2 CELESTIAL-FIXED AREA TARGETS AND SPACE 

VOLUMES 4-5 

5.0 DETAILED LOGIC FLOW 5-1 


APPENDIX A - SUBDIVIDING CONCAVE POLYGONS A-l 

APPENDIX B - CONIC INTERSECTIONS B-l 

REFERENCE R-l 



FIGURES 


Figure Page 

2-1 Aries mean-of-1950 coordinate system 2-2 

2-2 Rotating geocentric coordinate system 2-3 

2-3 Earth- referenced circle 2-5 

2-4 Geodetic coordinate system 

(a) Basic definitions 2-6 

(b) Detailed explanation 2-7 

2-5 Relationship between TEI and rotating geocentric 

systems 2-10 

2-6 Containment test for celestial circles 2-13 

2-7 Earth-referenced polygon 2-16 

2-8 Centroid of Earth- referenced polygons 2-18 

2-9 Containment test for Earth-referenced polygons . . 2-19 

2-10 Celestial polygon 2-22 

2-11 Containment test for celestial polygons 2-24 

2-12 Constant area polyhedron 2-26 

2-13 Boundary tests for constant area polyhedrons . . . 2-28 

2- 14 Containment test for constant area polyhedrons . . 2-30 

3- 1 Circumscribing polygons 3-6 

3-2 Circumscribing space volumes 3-8 

3-3 Intersection of celestial-fixed outer boundary 

cone with S/C orbit plane 3-12 

3-4 Earth-referenced outer boundary cone 3-15 

3-5 Spherical sector swept out by Earth-referenced 

outer boundary cone 3-16 

3-6 Definition of intersection points 3-18 

3- 7 Closest approach point 3-22 

4- 1 Functional flow for Earth-referenced area targets 

and space volumes 4-3 

4- 2 Functional flow for celestial-fixed area targets 

and space volumes 4-6 

5- 1 Detailed flowchart 5-2 


viii 



Page 


Figure 

A-l Examples of concave pentagons 

(a) One concave vertex A- 2 

(b) Two adjacent concave vertices A- 2 

(c) Two non adjacent concave vertices A- 2 

A- 2 Examples of subdividing concave pentagons 

(a) One concave vertex A- 3 

(b> Two adjacent concave vertices A- 3 

(c) Two nonadjacent concave vertices A- 3 

A- 3 Alternate subdivision A- 5 

B-l Intersection of two cones B-2 

B**2 Conditions when cones do not intersect ' 

(a) y > y , Y B “ 4 

1 T 2 

(b) y + smaller (Y t s Y ) < greater { y 1 : 

B-3 Spherical geometry to compute intersection 

points B-5 


ix 



ACRONYMS AND SYMBOLS 


AOS 

Acquisition-of-signal 

ATS VP 

Area Targets and Space Volumes Processor 

LOS 

Loss-of -signal 

MCC 

Mission Control Center 

M50 

Aries Mean-of-1950 

RNP 

Rotation, nutation, and precession 

S/C 

Spacecraft 

TCA 

Time of closest approach 

TEI 

True-of-epoch inertial 
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Semimajor axis 

"Adjustment factor" for computing altitude above 
the Fischer ellipsoid 

Vector along the centerline of a polyhedron 

Perpendicular distance to the side of a polyhedron 

Eccentricity 

Eccentric anomaly 

Flattening coefficient 

True anomaly 
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Inclination 

Mean anomaly 
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Unit outward normal vector 
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[RNP] 


TEI 
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U 
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♦ 
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y 

x 

y 


Radius of the Earth-referenced circle 
Mean equatorial radius 
Spacecraft position vector 

RNP matrix from the mean-of-1950 coordinate 
system to the true-of-epoch inertial system 

Time 

Epoch time corresponding to the RNP matrix 

Argument of latitude 

Right ascension 

Declination 

Geodetic latitude 

Geocentric latitude 

Central angle 

Longitude 

Earth gravitational constant 

Slant range vector 

Time of perigee passage 

Araument of perigee 

Earth rotation rate 

Right ascension of ascending node 

Ascending node vector 
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Subscripts: 

AOS 

B 

CA 

i 

LOS 

Superscripts: 

G 

TEI 

blank 


Denotes a parameter associated with the AOS point 

Denotes a vector measured from the lower boundary 
(i.e., base) of a polyhedron 

Denotes a parameter associated with the closest 
approach point 

Denotes the i t ^ 1 vertex or side of a polygon or 
polyhedron. Also denotes the i occurrence of 
an event 

Denotes a parameter associated with the LOS point 


Denotes a parameter in the rotating geocentric 
coordinate system 

Denotes a parameter in the true-of-epoch inertial 
coordinate system 

Denotes a parameter in the mean-of-1950 coordinate 
system 

Denotes a unit vector 
Denotes a vector 
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1.0 INTRODUCTION 


This document provides the level B/C mathematical specifications 
for the Area Targets and Space Volumes Processor (ATSVP) . Pur- 
suant to the requirements of reference 1, this processor is de- 
signed to compute the acquisition-of-signal (AOS) and loss-of- 
signal (LOS) times for the following: 

a. Area targets 

(1) Earth-referenced circles which are specified by a 
latitude, longitude, altitude, and radius. 

(2) Celestial circles which are specified by a right 
ascension, declination, and angular radius. 

(3) Earth-referenced polygons which are arbitrary Earth- 
fixed figures having up to five sides with the 
"corner points" defined by latitude, longitude, and 
altitude. 

(4) Celestial polygons which are arbitrary, inertially 
fixed figures having up to five sides with the corner 
points defined by right ascension and declination on 
the celestial sphere. 

b. Space volumes 

(1) Earth-referenced space volumes which are arbitrary. Earth- 
fixed polyhedrons having up to five sides. These volumes 
are defined by a lower-limit polygon at an altitude, h^, 
and the projection of this polygon to an altitude, hj. 

The corner points of the polygon are defined by 
latitudes and longitudes and rotate with the Earth. 

(2) Celestial-fixed space volumes which are arbitrary, inertial- 
ly fixed polyhedrons having up to five sides. These volumes 
are defined by a lower-limit polygon at an altitude, h^, 
and the projection of this polygon to an altitude, 

The corner points of the polygon are defined by right 
ascension and declination on the celestial sphere. 
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The AOS and LOS times for these targets are defined (ref. 1) as 
follows : 

a. Ground circles and polygons 

AOS - the time corresponding to the first subsatellite point 
to lie just inside the area. 

LOS - the time corresponding to the last subsatellite point 
just prior to exiting the area. 

b. Celestial circles and polygons 

AOS - the time corresponding to the first zenith point to lie 
just inside the area. 

LOS - the time corresponding tc the last zenith point just 
prior to exiting the area. 

c. Earth-referenced and celestial-fixed space volumes 

AOS - the time at which the spacecraft (S/C) is just entering 
the volume. 

LOS - the time just prior to the S/C exiting the volume. 

Six data tables will contain the information necessary to completely 
describe the area targets and space volumes. These tables (ref. 1) 
are as follows: 

a. Ground targets table containing 10 targets in 1 block of data. 

b. Celestial circles table containing 10 targets in 1 block of data. 

c. Ground polygons table containing 20 targets in 2 blocks of data. 

d. Celestial t olygons table containing 10 targets in 1 block of d?ta. 

e. Earth- referenced space volumes table containing 10 targets in 
1 block of data. 

f. Celestial-fixed volumes table containing 10 targets in 1 block 
of data. 

Section 2 of this document describes the characteristics of the 
area targets and space volumes and provides the mathematical 
equations necessary to determine whether the S/C lies within 
the area target or space volume. These equations provide a 
detailed model of the target geometry and will be used during 
the precise numerical search. 
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Section 3 discusses a sent: analytical technique for predicting 
the AOS and LOS time periods. This technique is designed to 
bound the actual visibility period using a simplified target 
geometry model and unperturbed orbital motion. Its principal 
purpose is to reduce the burden on the precise numerical search 
by eliminating regions of the S/C orbit where AOS and LOS times 
are physically impossible. Section 4 provides a functional over- 
view of the ATSVP. This section outlines the overall process 
required to determine precise AOS and LOS times. 

Section 5 presents the detailed logic flow for the ATSVP. This 
section integrates the functional overview presented in section 
4 with the equations and approach presented in sections 2 and 3 
and the appendixes. Appendix A discusses the procedure for sub- 
dividing complex concave polygons into two or more simpler convex 
segments. The purpose of this subdivision process is to pernit 
the equations in section 2 to be used on a segment-by- segment 
basis to test for containment. Appendix B provides a solution 
to the conic intersection equations used in section 3 for 
celestial-fixed targets. 


1-3 



2.0 AREA TARGETS AND SPACE VOLUMES CHARACTERISTICS AND CONTAINMENT 
CRITERIA 

The following subsections discuss the characteristics of each of 
the area targets and space volumes presented in section 1. The 
mathematical equations necessary to determine whether the S/C lies 
within the area target or space volume are also developed and dis- 
cussed. 

Two reference coordinate systems will be used. The inertial 
Aries mean-of-1950 (M50) coordinate system (fig. 2-1) will be the 
reference system when dealing with area targets and space volumes 
which remain inertially fixed. The rotating geocentric coordinate 
system (fig. 2-2) will be used when dealing with area targets and 
space volumes which rotate with the Earth. 

Each of the following six subsections is further subdivided into 
three topics: 

a. Procedure - a brief description of the steps to be performed. 

b. Equations - a statement of the input parameter requirements 
and development of the mathematical equations. 

c. Assumptions and limitations - a description of any simplifying 
assumptions and/or mathematical restrictions. 

For convenience, section 2.7 summarizes the equations for all of 
the area targets and space volumes. 
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I 


NAME; Aries mean-of-1950, Cartesian, coordinate system. 

ORIGIN: The center of the Earth. 

ORIENTATION: The epoch is the beginning of Besselian year 1950 or Julian 

ephemeris date 2433282.423357. 

The plane is the mean Earth's equator of epoch. 

The X^| axis is directed towards the mean vernal equinox of epoch. 

The Zj. axis is directed along the Earth's mean rotational axis 
of epoch and is positive north. 

The axis completes a right-handed system. 

CHARACTERISTICS: Inertial, right-handed, Cartesian system. 


Figure 2-1.- Aries mean-of-1950 coordinate system. 
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NAME: 

ORIGIN: 

ORIENTATION: 


CHARACTERISTICS: 


Geocentric coordinate system 
Center of the Earth 

X„ - Y„ plane is the Earth's true-of-date equator 
b b 

X G passes through the Greenwich meridian 
Z_ is along the Earth's rotational axis 
Y_ completes the right-handed system 

b 

Rotating, right-handed. Earth-fixed 


Figure 2-2.- Rotating geocentric coordinate system. 
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2.1 EARTH- REFERENCED CIRCLES 

Earth-referenced circles are defined to be circular ground target 
areas whose centers are defined by geodetic latitude, longitude, 
and altitude (fig. 2-3). The S/C lies within this ground target 
arc •. if its subsatellite point lies within the perimeter of the 
circular area. 

2.1.1 Procedure 

The following procedure will be used to determine whether the S/C 
lie- within the ground target area: 

a. The geodetic coordinates of the ground target area will be 
transformed to the geocentric system. 

b. The S/C position vector will be transformed from the M50 
system to the geocentric system. 

c. A test will be performed to determine whether the S/C subsat- 
ellite point lies within the perimeter of the ground target 
area. 

2.1.2 Equatio ns 

The following parameters are required: 

$ and X - geodetic latitude and longitude, respectively, of the 

center of the Earth- referenced circle (fig. 2-4) 

h - altitude of the Earth- referenced circle, measured with 

resoect co the Fischer ellipsoid of 1960 (fig. 2-4) 

r - arc .radius of the Earth-referenced circle (fig. 2-3) 
c 

R ^ - S/C position vector in the M50 system (fig. 2-3) 

SO 

t - time corresponding to R gc 
TEI 

r-SIP] - Rotation, nutation, and precession (RNP) matrix which 

is used to transform vectors from the M50 system to the 
true-of-epoch inertial (TEI) coordinate system 
t - epoch time corresponding to the RNP matrix 

R - mean equatorial radius for the Fischer elliosoid of 1960 

e.n 
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NAME: Geodetic coordinate system. 

ORIGIN: This system consists of a set of parameters rather than 

a coordinate system; therefore, no origin is specified. 

ORIENTATION: This system of parameters is based on an ellipsoidal model 

of the Earth (e.g. , the Fischer ellipse of 1960). For 
any point of interest we define a line, known as the geodetic 
local vertical, which is perpendicular to the ellipsoid 
and which contains the point of interest. 

h, geodetic altitude, is the distance from the point of 
interest to the reference ellipsoid, measured along the 
geodetic local vertical, and is positive for points out- 
side the ellipsoid. 

X is the longitude measured in the plane of the Earth's 
true equator from the prime (Greenwich) meridian to the 
local meridian, measured positive eastward. 

$ is the geodetic latitude, measured in the plane of the 
local meridian from the Earth's true equator to the geodetic 
local vertical, measured positive north from the equator. 

NOTE: A detailed explanation of declination, geodetic 

latitude, and geocentric latitude is provided 
on figure 2-4 (b). 

CHARACTERISTICS: Rotating polar coordinate parameters. Only position vectors 

are expressed in this coordinate system. Velocity vectors 
should be expressed in the Aries mean-of-1950 , or the Aries 
true-of-date, polar for inertial or quasi-inertial repre- 
sentations, respectively. The Fischer ellipsoid model 
should be used with this system. 

(a) Basic definitions. 

Figure 2-4.- Geodetic coordinate system. 
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Plane 

NAME: Geodetic coordinate system of point P. 

DEFINITIONS: h is the altitude of point P measured perpendicular 

from the surface of the referenced ellipsoid. 

4> is the geodetic latitude of point P. 

<|> c is the geocentric latitude of point P. 

6 is the angle between radius vector and equatorial plane 
(declination) . 

X is the longitude of point P. Angle (+ east) between 
plane of the figure and the plane formed by the Greenwich 
meridian. 


(b) Detailed explanation. 
Figure 2-4.- Concluded. 
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F - flattening coefficient for the Fischer ellipsoid of 1960 

w - Earth rotation rate 
e 


The first step is to transform the ground target area from the geo- 
detic system to the geocentric system. The vector from the 
center of the Earth to the center of the ground target area can 
be expressed in the rotating geocentric system by 


where 


B 


(h + a p ) cos X cos <j> 

I 

(h + a F ) sin X cos <J> 

[[h + (1 - F) 2 ap] sin <f>] 


a F 


em 


^:os 2 <J> + (1 - F) 2 sin 2 <j/ 


( 2 - 1 ) 


( 2 - 2 ) 


The second step is to transform the S/C position vector from the 
M50 system to the geocentric system. This is accomplished by 

a. Transforming the vector from the M50 system to the TEI system 
using the RNP matrix. 

b. Transforming the resultant vector from the TEI system to the 
geocentric system. 

The S/C position vector in the TEI system is given by 


S TEI - [bnp]™5 


sc 


sc 


(2-3) 
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The S/C position vector in the geocentric system (fig. 2-5) is 
given by 


R G = 
sc 


cos AX sin AX 0 
■sin AX cos AX 0 
0 0 1 


+TEI 

sc 

ITEI 


where 


AX = <o (t - t ) 
e e 


(2-4) 


(2-5) 


The final step is to determine whether the S/C subsatellite point 
lies within the ground target area (fig. 2-3). The angle from 
to the perimeter of the circular ground target area, y a » 
is given (in degrees) by 


Y A = 


m 


180 

TT 


The angle from to the S/C, y , is given by 


Y S 58 


= cos 


-1’ 


tf 3 

sc 

[Fi 

sc 


«S 


The S/C lies within the ground target area if 

Y S - y a 


(2-6) 


(2-7) 


(2-8) 
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Figure 2-5.- Relationship between TEI and 
rotating geocentric systems. 
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2.1.3 Assumptions and Limitations 

The following assumptions are implicit in the equations presented 

in section 2.1.2: 

a. The S/C geocentric subsatellite point is used to compute entry 
into the ground target area. 

b. The effects of polar nutation and precession from time t to 
time t can be neglected. 

c. The radius of the Earth- referenced circle, r , is a 
segment of arc. 
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2.2 CELESTIAL CIRCLES 

Celestial c.'rcles are defined to be circular areas on the celestial 
sphere (fig. 2-6). The center of this area target is defined by 
right ascension and declination. The criterion for a S/C to lie 
within this area is for the S/C zenith point to lie within the 
perimeter of the celestial circle. 

2.2.1 Procedure 

The following procedure will be used to determine whethe S/C 

lies within the celestial circle: 

a. The unit vector along the centerline of the celestial circle 
will be computed in the M50 system. 

b. The dot product between this vector and the S/C position 
vector will be formed to determine whether the S/C lies 
within the celertial circle. 


2.2.2 Equations 

The following parameters are required: 

a A and 6 A - the right ascension and declination, respectively, of 
the centerline of the celestial circle expressed in 
the M50 system (fig. 2-6). 

Y a - the celestial circle angular radius (fig. 2-6) . 

$ - S/C position vector in the M50 svstem. 

SC 


The unit vector from the center of the Earth to the center of the 

A 

celestial circle, C fi , is given by 


A 



cos a A cos 6^ 


sin a A cos 


sin o A / 


(2-9) 
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Centerline of 
celestial circle 



Figure 2-6.- Containment test for celestial circles. 


The angle between this vector and the S/C, Yg , is 


Y s * cos 


-1 




The S/C lies within the celestial circle if 


2.2.3 Assumptions and Limitations 


( 2 - 10 ) 


(2-li) 


None. 
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2.3 EARTH- REFERENCED POLYGONS 

The Earth- referenced polygon is defined to be an arbitrary planar 
figure having up to five sides (fig. 2-7) . This figure is fixed 
with respect to the rotating Earth. The corner points (i.e. , 
vertices) of this polygon are defined by geodetic latitude , longi- 
tude, and altitude. The basic criterion for penetration is to 
ensure that the S/C subsatellite point lies within the perimeter 
of the ground target area. 


2.3.1 Procedure 

The basic procedure for this ground target area is similar to the 

procedure presented in section 2.1.1. It consists of 

a. Transforming the geodetic coordinates of each polygon vertex 
to the geocentric system. 

b. Transforming the S/C position vector from the H50 system to 
the geocentric system. 

c. Testing to determine whether the S/C is interior to all 
planes defined by the sides of the polygon. 


2.3.2 Equations 

The following parameters are required: 


n 

and 



t 


[rnp] 


TEI 

M50 



number of sides 

geodetic latitude and longitude, respectively, of 
each vertex (fig. 2-7) 

geodetic altitude of each vertex (fig. 2-7) 

S/C position vector in the M50 system 
time corresponding to 5 

o C 

RNP matrix to transform from the M50 system to the 
TEI system 

epoch time corresponding to the RNP matrix 
mean equatorial radius for the Fischer ellipsoid of 1960 
flattening coefficient for the Fischer ellipsoid of 1960 
Earth rotation rate 
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Earth-referenced 

polygon 



Greenwich 

meridian 


Figure 2-7.- Earth-referenced polygon. 
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The first step is to transform the parameters defining each 
vertex of the polygon from the geodetic system to the rotating 
geocentric system. These vectors are given by 


where 



(h^ + ap) cos X^ cos 
(h^ + a F ) sin X^ cos 
JlK + (1 - F) 2 a p Jsin ^ 


i = 


1 » 2 , 3 , 


.n 


( 2 - 12 ) 


a_ is defined by equation 2-2 with replacing $. 

r X 

The S/C position vector in the geocentric system, ^g C » is then 
obtained by using equations 2-3 through 2-5. 

The next step in the procedure is to determine whether the S/C 
lies interior to all planes defined by the sides of the polygon. 
Figures 2-8 and 2-9 illustrate the geometry. All vectors in these 
figures are with respect to the geocentric system. The centroid 
of the ground target area, (fig. 2-8) is defined as 



*G _ i=l 
B n 


(2-13) 


The unit normal vectors to each side ot the polygon (fig. 2-9) are 
given by 


'G *1+1 * K 

1 l*Z +1 *"*ff 


i = 1,2,3 


-n-1 (2-14a) 



(2-14b) 


2-17 




Figure 2-8.- Centroid of Earth-referenced polygons 


2 -] 


Plane 2 



Center of 
the Earth 


Figure 2-9.- Containment test for Earth-referenced polygons. 
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The perpendicular distance from to each polygon side is 

d t = N? . i ^ 1,2, 3, . . .n (2-15) 

The S/C must lie interior to the i plane if 

N? . p® < d. (2-16) 

where 

- Kc - 5 B 

If equation 2-16 is satisfied for all sides, then the S/C sun- 
satellite point lies within the perimeter of the polygon. 

2.3.3 Assumptions and Limitations 

The following assumptions and limitations are implicit in the 
equations presented in section 2.3.2: 

a. The S/C geocentric subsatellite point is used to determine 
entry into the ground target area. 

b. The effects of polar nutation and precession from time 
to t can be neglected. 

c. The vertices of the polygon are specified in a counterclockwise 
order as viewed from the top. 

d. The polygon is convex; i.e. , the interior angles between the 
sides defining the vertices are less than 180 degrees.'*' 

e. The angular separation between two consecutive vertices is 
sufficient to permit a nonzero cross product (eq. 2-14). 

f. The S/C and the ground target area lie in the same hemisphere. 
The procedure discussed in section 3 will assure this 
condition. 

^Concave polygons can be accommodated by subdividing them into 
two or more convex polygons. Appendix A discusses this procedure. 
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2.4 CELESTIAL POLYGONS 

Celestial polygons are defined to be arbitrary figures having up 
to five sides with the corner points (i.e. r vertices) defined by 
right ascension and declination on the celestial sphere. The 
criterion for penetration into this area is to ensure that the 
S/C zenith point lies within the confines of the polygon. Figure 
2-10 illustrates the celestial polygon. It is noted that this 
polygon also represents a variable area polyhedron which remains 
inertially fixed. 

2.4.1 Procedure 

The basic procedure for this area target is similar to the proce- 
dure presented in section 2.3.1, i.e.,: 

a. The unit vector along the "centroid" of the polyhedron will 
be computed in the M50 system. 

b. Tests will be made to determine whether the S/C is interior 
to all planes defined by the sides of the polyhedron. 

2.4.2. Equations 

The following parameters are required: 
n - number of sides 

and 6^ - right ascension and declination, respectively, of each 

vertex in the M50 system (i = 1,2,3, ...n) 

$ - S/C position vector in the M50 system 

s c 


The first step is to determine the unit vector along the centroid 
of the polyhedron. The unit vectors from the center of the Earth 
to each vertex in the M50 system are given by 


A 



cos cos 6^ 
sin a. cos <5. 

l i 

sin 6. 

l 


x 1,2, 3, . . • n 


(2-18) 
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The unit vector along the centroid of the polyhedron is 


where 



( 2 - 19 ) 


(2-20) 


The final step is to determine whether the S/C lies interior to all 
planes defined by the sides of the polyhedron. Figure 2-11 illu- 

A 

strates the geometry. The slant range from C Lo the S/C is 



( 2 - 21 ) 


Equations 2-14 through 2-16 (with R. , C_, and p D replaced by 

X D D 

R., C, and p 0 , respectively) are then used to determine whether 

X 15 

the S/C is interior to all polyhedron planes. 


2.4.3 Assumptions and Limitations 

The following assumptions and limitations are implicit in the 

equations presented in section 2.4.2: 

a. The vertices of the polygon are specified in a counterclockwise 
order as viewed from the top. 

b. The polygon is convex; i.e., the interior angles between the 
sides defining the vertices are less than 180 degrees.'*' 

c. The angular separation between two consecutive vertices is 
sufficient to permit a nonzero cross product (eq. 2-14). 

d. The S/C and the area target lie within the same hemisphere. 

The procedure outlined in section 3 will assure this condition. 

Concave polygons can be accommodated by subdividing them 
into two or more convex polygons as discussed in appendix A. 
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2.5 EARTH- REFERENCED SPACE VOLUMES 

The Earth-referenced space volumes are arbitrary polyhedrons 
which are defined by a lower-limit polygon at an altitude, 
h^, and the projection of this polygon to an altitude, The 

vertices of this polygon are defined by geodetic latitude and 
longitude and rotate with the Earth. Figure 2-12 illustrates 
this type of space volume. As shown, the planar cross sectional 
area of this polyhedron remains constant with respect to altitude. 

2.5.1 Procedv e 

The following procedure will be used to determine whether the S/C 
lies within the space volume: 

a. The geodetic parameters defining each vertex of the lower 
boundary will be transformed to geocentric position vectors. 

b. The centroid vector to the lower boundary will be computed 
in the geocentric system. 

c. The S/C position vector will be transformed from the M50 
system to the geocentric system. 

d. Tests will be performed to ensure that the S/C lies above the 
lower boundary and belc : the upper boundary. If either of 
these tests fails, then no further computations are re- 
quired. 

e. Assuming the previous step is passed, tests will be performed 
to determine whether the S/C is interior to all planes 
defined by the sides of the polyhedron. 

2.5.2 Equations 

The following parameters are required: 
n - number of sides 

and A i - geodetic latitude and longitude, respectively, of 
each vertex of the lower boundary 
h^ - geodetic altitude of the lower boundary 

h 2 - geodetic altitude of the upper boundary 
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Space volume centerline 



Figure 2-12.- Constant area polyhedron. 
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sc 


[rnp] 


TEI 

M50 


em 


OJ 


e 


S/C position vector in the M50 system 
time corresponding to $ 

sc 

RNP matrix to transform from the M50 to TEI system 
epoch time corresponding to the RNP matrix 
mean equatorial radius for the Fischer ellipsoid 
of 1960 

flattening coefficient for the Fischer ellipsoid 
of 1960 

Earth rotation rate 


The first step is to transform the parameters defining each 
vertex of the lower boundary to geocentric position vectors. 
Equation 2-12 (with h^ = 0) provides the necessary transformation 


The centroid of these vectors is given by 



( 2 - 22 ) 


The geocentric vector to the centroid of tho lower boundary 
(fig. 2-13) is given by 


3 


C G + h. 



(2-23) 


Next, the S/C position vector in the geocentric system, , is 

computed using equations 2-3 through 2-5. 

The fourth step in the procedure is to ensure that the S/C lies 

between the upper and lower boundaries. Figure 2-13 illustrates 

the geometry. All vectors in this figure are with respect to 

the geocentric system. The slant range from c G to the S/C, p G 

D B 


2-27 




Figure 2-13.- Boundary tests for constant area polyhedrons. 
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is given by 


n - Ko - «*-*« 

The S/C lies above the lower boundary and below the upper boundary 
if 


h 2 " h l - 



(2-25) 


Assuming that equation 2-25 is satisfied, the final step is to 
determine whether the S/C lies interior to the side planes of 
the polyhedron. Figure 2-14 illustrates the geometry. The unit 
normal vectors from C_ to each side of the polygon are given by 

O 





1,2,3. . . n— 1 

(2-26a) 


(2-26b) 


Equations 2-15 through 2-17 are then used to determine whether the 
S/C is contained within the polyhedron. 


When dealing with concave space volumes, using the segmentation 
procedure discussed in appendix A, the following procedure must 
be follevod to prevent geometric distortion of the segments. 

•+G 

a. The centroid of the entire space volume is used for xn 

equations 2 24 and 2-25 to test the upper and lower boundary 
constraints. 
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b. The centroid of the entire space volume is also used for C? 
in equation 2-26 to construct the unit normal vectors on a 
segment- by- segment basis. 

c. The segment centroid is used in equations 2-15 through 2-17 
to test for parameter containment. 


2.5.3 Assumptions and Limitations 

The following assumptions and limitations are implicit in the 

equations presented in section 2.5.2: 

a. The altitude of the lower boundary, h^, is measured with 
respect to the centroid of the vertices which lie on the 
Fischer ellipsoid. 

b. The effects of polar nutation and precession from time t 
to t can be neglected. 

c. The vertices are specified in a counterclockwise order as 
viewed from the top. 

d. The planar area of the polyhedron is convex; i.e., the inter- 
ior angles between the sides defining the vertices are less 
than 180 degrees. 1 

e. The local horizon for determining whether the S/C is between 
the upper and lower boundaries is perpendicular to the centroid 
vector. 

f. The distance between two consecutive vertices is sufficient 
to permit a nonzero cross product (eq. 2-26) . 


Concave polyhedrons can be accommodated by subdividing them 
into two or more convex portions as discussed in appendix A. 
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2.6 CELESTIAL-FIXED SPACE VOLUMES 

The celestial-fixed space volumes are arbitrary polyhedrons 
which are defined by a lower-limit polygon at an altitude , 
h^, and a projection of this polygon to an altitude, l^- The 
vertices of the polygon are defined by right ascension and de- 
clination and are inertially fixed. Figure 2-12 can also be 
used to illustrate this type of space volume. As mentioned 
previously, the planar cross sectional area of this polyhedron 
remains constant with respect to altitude. 

2.6.1 Procedure 

The basic procedure for this space volume is very similar to the 
procedure presented in section 2.5.1. It consists of: 

a. Computing the centroid vector to the lower boundary in the 
M50 system. 

b. Testing to ensure that the S/C lies above the lower boundary 
and below the upper boundary. If either of these tests fails, 
then no further computations are required. 

c. Assuming the previous step is passed, further tests will be 
performed to determine whether the S/C is interior to all 
planes defined by the sides of the polyhedron. 


2.6.2 Equations 

The following parameters are required: 


n 

a . and 6 . 


em 


number of sides 

right ascension and declination, respectively, of 
each vertex of the lower boundary in the M50 system 
altitude of the lower boundary 
altitude of the upper boundary 
mean equatorial radius 
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The unit vector to the centroid of the polyhedron, C, is com- 
puted via equations 2-18 through 2-20. The vector to the centroid 
of the lower boundary is given by 

- ( R ea + h l) a « 2 - 27 ' 

-►Q 4Q 

Equations 2-24 and 2-25 (with variables Rr and replaced by 

-v sc B 

R and C_, respectively) are then used to ensure that the S/C 

SO D 

lies between the upper and lower boundaries. 

Assuming equation 2-25 is satisfied, equation 2-26 (with C fi 
and £? replaced by and R. , respectively) is used to 

define the unit normal vectors to each side of the polyhedron. 
Finally, equations 2-15 through 2-17 are used to determine whether 
the S/C lies within the space volume. 


When dealing with concave space volumes, using the segmentation 
procedure discussed in appendix A, the following procedure must 
be followed to prevent geometric distortions of the segments. 


a. The centroid of the entire space volume is used for C„ in 

O 

equations 2-24 and 2-25 to test the upper and lower boundary 
constraints . 

b. The centroid of the entire space volume is also used for C G 

B 

equation 2-26 to construct the unit normal vectors on a 
segment-by- segment basis. 

c. The segment centroid is used in equations 2-15 through 2-17 
to test for perimeter containment. 


in 
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2.6.3 Assumptions and Limitations 

The following assumptions and limitations are implicit in the 

equations presented in section 2.6.2: 

a. The altitudes of the lower and upper boundaries are measured 
with respect to the mean equatorial radius. 

b. The vertices are specified in a counterclockwise order as 
viewed from the top. 

c. The planar area of the polyhedron is convex; i.e. , the in- 
terior angles between the sides defining the vertices are 
less than 180 degrees.'*' 

d. The distance between two consecutive vertices is sufficient 
to permit a nonzero cross product (eq. 2-26) . 


Concave polyhedrons can be accommodated by subdividing 
them into two or more convex oortions as discussed in appendix A. 
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2.7 SUMMARY OP EQUATIONS 

This section summarizes all of the equations presented in the 
previous sections for the various area targets and space volumes. 
The order of presentation and equation numbers correspond to the 
computation sequence discussed in the text. 


2.7.1 Earth- Referenced Circles 


where 



(h + a cos A cos 4> 

• (h + a^j sin A cos $ 
[h + (1 - F) 2 a plsin 

k * j 


( 2 - 1 ) 


a F 35 


em 


y, 2, . ,, ZT2 _T_27 


( 2 - 2 ) 


' cos 4> + (1 - F) sin (j> 


t>TEI s FRNPl™ ^ 

sc LRNPJ M50 R sc 


(2-3) 



cos A A 
-sin AA 
0 


sin AA 
co AA 
0 


0 

0 

1 


±TEI 

R sc 


t’TEI 


(2-4) 


where 



(2-5) 
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Ys i 


(2-6) 

( 2 - 7 ) 

(2-8) 
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where 


where 


■p 


em 




( 2 - 2 ) 


R TEI = rRNPl TEI £ 
R sc l «np j m5o r sc 


(2-3) 


R sc 


cos AA sin AA 0 
•sin AA cos AA 0 
0 Oil 


+TEI 

R sc 

|TEI 


( 2 - 4 ) 


AA 


■ “e (* - *.) 


( 2 - 5 ) 


±n 


? G = i=1 
B n 


(2-13) 


go tf+l x S 
1 = lS? +1 x *?l 


i = 1,2,3, — n-1 (2-14a) 


N = 
n 


d. 

l 


8? x 8 s 

1 n 


IS? x i£i 

»? • (*? - S) 1 


(2-14b) 


= 1,2,3,. . .n (2-15) 


»? • Pb i a i 


i = 1,2,3, . . .n (2-16) 


2-37 



where 


*2 ■ - *S 


(2-17) 


2.7 .4 Celestial Polygons 


where 


R i = 


cos a . cos . 
i x 

sin a. cos 6. 
i x 

sin 6^ 


i = 1,2,3, ...n (2-18) 


C * 


1*1 


(2-19) 


n 


E a 

R i 

2 = 1=1 


n 


( 2 - 20 ) 


- he - £ 


( 2 - 21 ) 


N. = 


R., X R. 
X+l X 

|R i+ l X R.| 


i = 1,2,3, . . .n-1 (2-14a) 


N = 
n 


R. X R 
1 n 

a ~ i 

R, X R 
1 n 1 


(2-14b) 


! i = ’ ( R i ' i = 1,2,3, ...n (2-15) 


Si • P B < d 4 


i * 1,2,3, ...n (2-16) 
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2.7.5 Earth-Referenced Space Volumes 


where 


a F cos X^ cos 4>^ 


*G 


~ « a F sin X^ cos 4 k 

(1 - F) 2 a p sin 4> i 


i * 1,2,3, 


a_ - 
F 


em 




E *? 


$ G a* 


n 


s' 3 
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+TEI 

sc 

= 

[rnp] tei 

lkw ^ j M 50 

sc 


cos 

AX 

sin AX 

O' 

s 

-sin 

AX 

cos AX 

0 


0 


0 
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IG 

j^TEI 
SC 

fTEI 


where 


AX 


" “e (‘ - S) 


P° = R° - 2S 

B SC B 


( 2 - 12 ) 

(2-2) 

( 2 - 22 ) 

( 2 - 23 ) 

( 2 - 3 ) 

( 2 - 4 ) 

( 2 - 5 ) 

( 2 - 24 ) 


2-39 



h 2 - h, > p 


-►G 


g B * ( ' - g i,l> 


l«B X 


r g 

i+I 


g G g * (*n - tf ) 
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2.7.6 Celestial-Fixed Space Volumes 


cos cos 6^ 
R. = . sin a . cos 6 . 

l li 

sin 6 . 
i 


■ ^ 

E R i 




C = (r + h.) C 
B v em 1/ 


2-40 



(2-24) 




(2-26b) 

= 1,2 ,3, . . .n (2-15, 

= 1,2,3, .. .n (2-16) 
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. 0 SEMIANALYTICAL ALGORITHM TO PREDICT AOS AND LOS TIMES 


This section presents a semianalytic technique for predicting AOS 
and LOS times. This technique is designed to bound the actual 
visibility periods and uses a simplified geometry model in con- 
junction with unperturbed orbital motion* Its principal purpose 
is to reduce the burden on the precise numerical search by 
eliminating r^" Lons of the S/C orbit where AOS and LOS times are 
physically impossible. 

The basic procedure consists of the following: 

a. Determining the parameters of an outer boundary cone which 
circumscribes the entire area target or space volume. 

b. Determining whether the orientation of the S/C orbic plane 
will permit any intersections of the outer boundary cone 
with the S/C orbit plane. 

c. Assuming intersections are possible for the following: 

(1) Celestial-fixed targets 

(a) Compute the conic intersection points of the outer 
boundary cone with the S/C orbit plane. 

(b) Compute predictions of AOS and LOS times 

( fc AOS and fc LOS ' respectively) for each 
K K 

revolution based upon these intersection points. 

(2) Earth-fixed targets 

(a) Compute the time of closest approach (TCA) of the 
S/C position vector with the centerline of the 
outer boundary cone and the corresponding closest 
approach angle (y ) . These computations are 
performed on a revolution-by-revolution Lisis. 

(b) If y ca is less than or equal to the half-cone 
angle of the outer boundary cone, then compute 
predictions of AOS and LOS times based upon the 
times when the S/C lies on the perimeter of the 
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outer boundary cone. 


’AOS. 


and t 


LOS. 


These times are denoted by 


After the AOS and LOS time predictions are determined, the precise 
numerical search (using the detailed area target and space 
volume geometry model discussed in sec. 2) can then be performed 
over the reduced time regions from t._ e to t__ e . 


The following parameters are required to predict the AOS and LOS times: 


fc start 

fc end 

At 



V 

sc o 


[RNP] 


TEI 

M50 


t 

e 


- start time for the search 

- end time for the search 

- time increment for the search. For the purposes of 
this algorithm. At will be used as a tolerance for 
converging upon the AOS and LOS times for Earth- 
referenced targets. 

- S/C position vector in the M50 system at t start * 

- S/C velocity vector in the M50 system at t 

S LaLL 

- RNP matrix to transform from the M50 to the TEI 
system. 

- epoch time corresponding to the RNP matrix. 


The following subsections provide the necessary equations. 

Section 3.1 discusses the computation of the outer boundary cone 
parameters. Section 3.2 discusses the technique to determine 
whether any AOS or LOS times are possible. Section 3.3 presents 
the equations to predict AOS and LOS times. Section 3.4 discusses 
the assumptions and limitations implicit in these equations. 

Two reference coordinate systems are used in these computations. 
The M50 coordinate system is used when dealing with celestial- 
fixed targets and the TEI system is used for Earth-referenced 
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targets. The TEI system was selected for Earth-referenced targets 
in order to model the effects of polar nutation and precession 
on the Earth's spin axis. The 2 axis of the TEI system is aligned 
with the Earth's spin axis at the epoch time of the BNP matrix. 

Tor celestial-fixed targets, the M50 Keplerian elements (a, e, 

i, ft, w and M ) at the start time will be computed from R 

o •'*o 

and V 8C * 


For Earth-referenced targets, the TEI Keplerian elements at 

4. irgT "^TEI • 

start time will be computed from R and V where 

sc 0 s 0 

the 

- C ra,p £ I S sc 0 


(3-1) 

? sc 0 

0 0 


(3-2) 
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3.1 COMPUTING OUTER BOUNDARY CONE PARAMETERS 


The outer boundary cone is defined by 

a - right ascension of the centerline 
o 

d> - latitude of the centerline 
o 

Y a - half-cone angle 


For Earth- referenced targets these parameters will be defined in 
the TEI system at ^gtart*^ For ce l est i a l~fi xe< ^ targets, these 

parameters will be defined in the M50 system and remain invariant 
with respect to time. 


3.1.1 Earth-Referenced Circles 

The unit vector along the centerline of the cone is given in the 
TEI system by 


C 


TEI 


cos AX 
sin AX 
0 


sin AX 
cos AX 
0 


0 

0 

1 


TEI 


% 


I I 


(3-3) 


where 

C_ is given by equation 2-1 
B 

AX is given by equation 2-5 with t equal to t . . 

s t3ru 


Thus, the right ascension and latitude of the centerline is given 
in the TEI system (at the start of the search) by 


tan 


-1 


r TEI 

5l_ 

C TEI 

x 


(3-4) 


The TEI system is selected for Earth-referenced targets in 
order to model the polar nutation and precession from M50 reference 
to the epoch time of the RNP matrix. 


3-4 



(3-5) 


. . -1 / _TEll 

« sin 1C I 
o ' g ' 

where 

"TEI ^ Lt BI _TEI _TEI 

\ X ' c y ' c z 

The half-cone angle, y A , is given by equation 2-6. 

3.1.2 Celestial Circles 

The right ascension, latitude, and half-cone angle in the H50 
system arc- given via input 


“o - “a 

(3-6) 

♦o*« A 

(3-7) 


3.1.3 Earth- Referenced Polygons 

The unit vector along the centerline of the outer boundary cone is 
given in the TEI system by equation 3-3 where equation 2-13 is used 
to compute C^. Equations 3-4 and 3-5 are then used to compute 
a Q and <j> Q , respectively. The outer boundary cone is centered 
along the centroid vector and has an angular radius which contains 
the vertex furthest from the centroid vector (fig. 3-1). Thus, 
the half-cone angle is given by 


where 


max 



r9 . 

cS 


cos * 

1 

B 


|r9| 

i*Ji 


(3-8) 


max implies the maximum algebraic value 
q is given by equation 2-12. 
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Figure 3-1.- Circumscribing polygons. 


3.1.4 Celestial Polygons 


The unit vector along the centerline of the celestial polygon is 

^tei 

given by equation 2-19. Equations 3-4 and 3-5 (with C replaced 
by C) are then used to compute <x Q and <J> Q , respectively. The 
half-cone angle which contains the vertex furthest from the centroid 
is given by 


Y a = max 


cos | • C | ^ lf2|3f*. . nj 


(3-9) 


where 


R. 

i 


is given by equation 2-18. 


3.1.5 Earth-Referenced Space Volumes 

The unit vector along the centerline of the outer boundary cone 
for Earth- referenced space volumes is given by equation 3-3 with 

6 G replaced by $ G (eq. 2-22). Equations 3-4 and 3-5 
B 


are then used to compute 


<* o and <f> Q , respectively. Since the 


Earth- referenced space volume represents c constant area polyhedron 
(sec. 2.5), the outer boundary cone to be used for predicting AOS 
and LOS times circumscribes the space volume planar area at the 
S/C altitude. However, for noncircular orbits, the S/C altitude 
varies as a function of time and thus the size of the outer 
boundary cone would also vary as a function of time. The size of 
the cone is largest at perigee and smallest at apogee. Since the 
purpose of the outer boundary cone is to bound the actual visibility 
region, the half-cone angle at the S/C perigee altitude will be 
computed and held constant with respect to time. The vector from 
the center of the Earth to each vertex of the space volume pro- 
jected to the S/C perigee altitude, R^ , (fig. 3-2) is given by 


= a? 


« P - i«r 


C B 


I * 


B 


i = 1,2,3, . . .n (3-10) 
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Figure 3-2.- Circumscribing space volumes. 
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where 


Q 

K? is given by equation 2rl2 
C fi is given by equation 2-23 


"p 


is the S/C perigee radius 
■ a (1 - e) 


(3-11) 


Thus the half-cone angle which contains the vertex furthest from 
the centerline of the cone is given by 


Y a = max 


cos 


-1 


n' 


*B 


l*fl I^b I 


i ™ 1,2,3, .. .n 


(3-12) 


3.1.6 Celestial-Fixed Space Volumes 

The unit vector along the centerline of the outer boundary cone 
for celestial- fixed space volumes is given by equation 2-19. 

^tei * 

Equations 3-4 and 3-5 (with C replaced by C) are then used 
to compute a Q and <t» 0 , respectively. Since the celestial- fixed 

space volume also represents a constant area polyhedron (sec. 2.6), 

the outer boundary cone to be used for predicting AOS and LOS times 

circumscribes the space volume planar area at the S/C perigee 

altitude (sec. 3.1.5). The vector from the center of the Earth to 

each vertex of the space volume at the S/C perigee altitude , 

* . 

R i , given by 



R. + 
l 




(3-13) 
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where 


Thus 


A 

R^ is given by equation 2-18 
6.. is given by equation 2-27 

Jd 

R is the S/C perigee radius 
P 

(equation 3-11) 



(3-14) 
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3.2 DETERMINING WHETHER AOS AND LOS TIMES ARE POSSIBLE 


The following computations will determine whether the relative 
orientation of the S/C orbit plane and the outer boundary cone 
will permit any intersections between these two figures. If no 
intersections are possible, then no AOS or LOS times are possible. 
Since the Earth- referenced targets rotate with the Earth and the 
celestial-fixed targets remain inertially fixed, two different 
types of tests are required. The celestial-fixed situation is 
the simplest and is discussed first. 

3.2.1 Celestial-Fixed Targets 

Figure 3-3 illustrates a celestial-fixed outer boundary cone and 
a S/C orbit plane. The S/C orbit plane is perpendicular to the 
unit angular momentum vector which is given by 


h = 



(3-15) 


where 

are the S/C position and velocity vectors, re- 
spectively, in the M50 system at the beginning 
of the search. 

A 

The centerline of the outer boundary cone, C, is given in the 
M50 system by 


sc. 


and V 


sc. 



cos 

a cos 

o 

*o 

35 < 

sin 

a cos 

o 

♦o 



sin <f> 

o 

1 


(3-16) 
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Figure 3-3.- Intersection of celestial-fixed 
outer boundary cone with S/C orbit plane. 
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where 

a Q and <J> Q were computed using the procedure discussed in 
section 3.1. 

The intersection points of the outer boundary cone with the S/C 

A /V 

orbit plane are denoted by 1^ and on figure 3-3. These 

vectors can be found by solving the following conic intersection 
equations* 

h • I = 0 (3-17a) 

A A 

C • I * cos y a (3-i7b) 

1 - 1=1 (3-17c) 

A 

These equations provide 0 f 1. or 2 solutions for I. If no 
solutions are produced, then the relative orientation of the 
outer boundary cone and the S/C orbit plane does not permit 
intersections, and no AOS or LOS times are possible. If only 
one solution is produced, the AOS and LOS times are the same. 

If two solutions are produced, the AOS and LOS times predictions 
are made using the method outlined in section 3. 3. 


The method for solving these equation- is presented in 


appendix B. 
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3.2.2 Earth-Referenced Targets 

Figure 3-4 illustrates an Earth- referenced outer boundary cone 
at a particular instant of time. This cone rotates about the 
Earth's polar axis. Hence, over one sidereal day. the volume of 
space swept out by the outer boundary cone describes a spherical 
sector whose upper and lower limits are determined by y a and 
<J> o (fig. 3-5). Intersections of the S/C orbit plane with this 
volume of space are possible only if the S/C orbit inclination 
exceeds the lower limit of the spherical sector. 

For posigrade and polar orbits, the S/C orbit plane will intersect 
the spherical sector only if^ 

i>|0 o |-Y A ( 3-18a) 

For retrograde orbits, intersections are possible only if ^ 

180 - i > U 0 | - Y a ( 3—1 8b) 


where 

i is the S/C orbit inclination in the TEI system 

If equation 3-18 is not satisfied, then no intersections are 
possible and thus no AOS or LOS times are possible. If equation 
3-18 is satisfied, then the approach presented in section 3.3 is 
used to predict AOS and LOS times. 


^The absolute value sign is used in these equations to ac- 
commodate targets in both the northern and southern hemispheres. 
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Figure 3-4.- Earth-referenced outer boundary cone. 
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3.3 PREDICTING AOS AND LOS TIMES 

The following subsections present the equations necessary to 
predict AOS and LOS times for celestial-fixed and Earth-referenced 
targets. These equations assume that the tests performed in 
section 3.2 have indicated that intersections of the outer boundary 
cone with the S/C orbit plane are possible. 


3.3.1 Celestial-Fixed Targets 

A A 

The intersection vectors, 1^ and I (sec. 3.2.1), can be 
used directly to predict the AOS and LTS times. Figure 3-6 
illustrates the geometry. The unit vector in the direction of 

A 

the ascending node, ft, is given by 


ft 


cos ft 
sin ft 
0 


(3-1S) 


where 

ft = right ascension of the ascending node in the M50 
system 

A 

The angles between ft and the two intersection vectors are given 
by 


u 


1 


tan 


A A A 


h • (« X 1^ 



u 


2 


tan 


h • (ft X I 2 ) 

A. A 

ft • I 2 


(3-20a) 


(3-20b) 
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One of these angles corresponds to the AOS point and the other to 
the LOS point. The appropriate angle can be determined by noting 
that the visibility period must be less than or equal to one half 
of the orbit. 


Thus let 


If K > 0 


If K < 0 


K — tan” 1 


h . (I 1 X I 2 ) 


/V A 



(3-21) 


U A0S 

U 1 

(3-22a) 

U LOS 

u 2 

(3-23a) 

U AOS 

U 2 

(3-22b) 

U L0S 

U 1 

( 3— 23b) 


The true anomalies of the AOS and LOS points are given by 


f AOS/LOS U AOS/LOS * “ 


(3-24) 


where 

w is the argument of perigee of the S/C orbit 
These can be converted to AOS and LOS times by 


* _ M A0S/L0S , 

AOS/LOS 

M 


(3-25) 


where 


M A0S/L0S ” E AOS/LOS “ e sin E A0S/L0S (3-26) 
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e aos/los 


(3-27) 


-»-- l |rfRT *» ( f ^)| 


(3-28) 

(3-29) 
(3-30) 

R is the period number measured with respect to 
fc start 

a is the semimajor axis of the S/C orbit 

e is the eccentricity of the S/C orbit 

M is the S/C mean anomaly at t . 
o start 

y is the gravitational constant 


M = 


360 

-^r- (in degrees per unit of time) 


■J? 


T - Start - r 2 + < R * » T 

M 


Since both the orbit plane and the celestial-fixed target remain 
inertially fixed, u ^q S and u los rema i n constant. Thus» the 
AOS and LOS times for future revolutions differ by only the 
orbital period. The AOS and LOS times for the subsequent N 
revolutions are given by 


t AOS/LOS K+1 ~ fc AOS/LOS K + T 


K = 1,2,3. ..N (3-31) 
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3.3.2 Earth-Referenced Targets 

Computation of AOS and LOS times for Earth-referenced targets is 
somewhat more complicated than for celestial- fixed targets since 
the orientation of the Earth-referenced target varies as a func- 
tion of time. The following procedure will be used to estimate 
the AOS and LOS times: 

a. The TCA of the center of the Earth-referenced outer boundary 
cone with the S/C position vector will be determined along 
with the corresponding closest approach angle# Yca . These 
computations will be performed on a revolution-by-revolution 
basis. 

b. The closest approach angle from the previous step will be 
compared with y R . If y ca > Y , then no AOS or LOS times 
are possible for that revolution, and the previous step will 
be repeated for the next revolution. If Yca £ Y A » then 
AOS and LOS times will be computed by solving for the times 
when the angle between the center of the outer boundary cone 
and the S/C position vector is equal to y^. 

3. 3. 2.1 Determining TCA 

The closest approach point occurs when the S/C position vector 
and the center of the outer boundary cone lie in the same hemis- 
phere and are coplanar (fig. 3-7) . Both the S/C and the center 
of the outer boundary cone are moving with time. The center of 
the cone is rotating about the Earth's Z axis at the Earth's 
rotation rate. The S/C is moving about the orbit angular momen- 
tum vector at a rate which is dependent upon both the orbit period 
and eccentricity. As a result, the TCA cannot be analytically 
determined for noncircular/nonequatorial orbits. The TCA can, 
however, be determined iteratively by the following: 
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Figure 3-7.- Closest aooroach print. 
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a. Estimating the closest approach point (for the first revolu- 
tion, this estimate would be based upon t . . ; for sub- 

start 

sequent revolutions, this estimate would be based upon the 
TCA of the previous revolution plus one orbital period) . 

b. Determining the time at which the S/C will reach this point. 

c. Updating the closest approach point based upon the time 
estimate from the previous step. 

d. Repeating the two previous steps until the difference in 
time between two successive iterations is less than or equal 
to the time step. At. 


The following discussion develops the equations necessary to 
compute the TCA. 


The unit vector from the center of the Earth to the center of 
the outer boundary cone is given by 


cos a cos <|> 

o 

sin a cos <J> (3-32) 

o 

sin d> 
r o 


where 



a 


+ “e (t - ‘start 1 


(3-33) 


a is the right ascension of the centerline at t . , 

o start 

(sec. 3.1) 


“> e is the Earth rotation rate 

<|> o is the latitude of the centerline (sec. 3.1) 
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The TCA, is then given by 

f CA = U CA ” a ' 



( 3 - 36 ) 

( 3 - 37 ) 

( 3 - 38 ) 

( 3 - 39 ) 
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where 


t is given by equation 3-30 
• 

M is given by equation 3-28 
e is the orbit eccentricity 
w is the argument of perigee 

This time would then be used in equations 3-32 and 3-33 to update 
the value of C. Equations 3-35 through 3-39 would be repeated 
until the value of t CA between two successive iterations is 
less than or equal to the time step. At. After the TCA has been 
found, the corresponding closest approach angle, Y CA * would be 
computed. 


The unit vector from the center of the Earth to the S/C is given 
by 


r = 


cos Q cos u - sin Q sin u cos i 
sin cos u + cos Q sin u cos i 
sin u sin i 


(3-40) 


The closest approach angle is given by 

-1 


Yca = cos 


r CA * C CA 


where 


(3-41) 


r CA is 9 iven by equation 3-40 evaluated at u_. 
C CA is 9iven by equations 3-32 and 3-33 evaluated 

at t CA* 
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If y ca is greater than the outer boundary cone angle, y a 
( sec. 3.1), then no AOS or LOS times are possible for this 
revolution. In this event, the procedure outlined is repeated 
for the next revolution. The initial estimate of TCA for the 
next revolution would be equal to the TCA for the current re- 
volution plus one orbital period. The process continues until 
either £ Y n or the TCA exceeds the end time of the search, 

Ln, A 


'end' 


If 


Y^ a 5 y a » estimates of AOS and LOS times are made 


using the method outlined in section 3. 3. 2. 2. 


3. 3.2.2 Determining AOS and LOS Times 

After a feasible' TCA has been found, the AOS and LOS times can 
be computed by solving the following equations. 

A A 

r • C = cos y a (3-42) 

Substituting equations 3-40 and 3-32 into this expression and 
simplifying yields 


cos (ft - a) cos <{> cos u - sin (ft - a) cos <J> Q cos i sin u 
+ sin <p sin i sin u = cos y. (3-43) 

O H 


Equation 3-43 cannot be explicitly solved for time for noncircular/ 
nonequat.orial orbits. However, this equation can be solved 
numerically using the method of successive substitution. For 
convenience, equation 3-43 can be rewritten as 


A n cos u + B 1 sin u = ^ 

where 


> 

II 

COS 

$ 

y o 

cos 

(ft 

- a) 

B l " 

sin 

T o 

sin 

i - 

cos d> cos i sin (fi - a) 
o 

c l = 

cos 

y a 





(3-44) 


(3-45) 

(3-46) 

(3-47) 
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Assuming a is constant between iterations, equation 3-44 has 
the following solution 


u AOS/LOS 


_1 

tan 




(3-48) 


The minus sign will produce AOS solutions and the plus sign will 
produce LOS solutions. The procedure is as follows: 


a. Solving for AOS times, fc A os 

(1) Use initial estimate of t AQS = t CA . 

(2) Compute A B and Cj^ using equations 3-33, 3-45, 
3-46, and 3-47. 

(3) Compute u c using equation 3-48 with the minus sign. 

(4) Compute revised estimate of AOS time using equations 
3-24 through 3-30. 

(5) If the change in time from the previous iteration is 
less than At, then a solution has been found. Other- 
wise, repeat steps (2) through (5) with the revised time 
estimate. 

b. Solving for LOS time, t LQS 

The solution for LOS time is similar to the procedure for 

AOS time except that 


(1) The initial estimate of LOS time in step a(l) is given 
by 


( 2 ) 

(3) 


fc LOS = fc CA 
u^Qg is computed in 
with the plus sign. 
The revised estimate 


+ (t CA " ^OS 5 (3-49) 

step a (3) by using equation 3-48 

of LOS time is produced in step a (4). 
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The computations discussed in this section are performed for 
each revolution in which was found to be less than or equal 

to y a - In the special case where y ca is exactly equal to y a » 
the AOS and LOS times are identical and equal to the TCA. Hence 
computation of AOS and LOS times c?n be bypassed. 
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3.4 ASSUMPTIONS AND LIMITATIONS 

The following assumptions are implicit in the equations presented 
in sections 3.1 through 3.3. 

a. Unperturbed orbital motion over the time period t ^ i t 
b start. 

b. The effects of polar nutation and precession from t to 

fc end can be ne< ? lecte< 3* 
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4.0 FUNCTIONAL OVERVIEW OF AOS AND LOS TIME COMPUTATIONS 

This section provides a functional overview of the ATSVP. The 
overax 1 procedure for determining precise AOS and LOS times is 
presented and discussed. The basic procedure consists of 
predicting the AOS and LOS times using the algorithm presented 
in section 3 and then refining these predictions by performing 
a sequential time search using the equations presented in section 2. 

The AOS times are defined to be the time point at which the S/C 
first enters the area target or space volume. If the S/C lies 
within the area target or space volume at the beginning of the 
user-specified search period (t s t a rt^ * the f ^ rst A0S time will 
be set equal to the start time. LOS times are defined to be 
the last time point prior to the S/C exiting the area target 
or space volume. If the S/C lies within the area target or space 
volume at the end of the user-specilied search period (t , then 
the last LOS time will be set equal to the end time. It should 
be noted that the time increment used to perform the sequential 
time search will limit the accuracy and resolution of the AOS 
and LOS times (e.g., if the time increment is 1 minute, this 
implies that the AOS and LOS times will be determined to the 
nearest minute and that visibility periods of less than 1 minute 
may be skipped) . 

The area targets and space volumes defined in section 1 fall into 
two general categories: 

a. Earth-referenced which include 

(1) Earth-referenced circles. 

(2) Earth-referenced polygons. 

(3) Earth-referenced space volumes. 

b. Celestial-fixed which include 

(1) celestial-fixed circles. 

(2) celestial-fixed polygons. 

(j) celestial-fixed space volumes. 
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The logic for computing the AOS and LOS times for each of these 
categories is presented in the following subsections . 


4.1 EARTH- REFERENCED AREA TARGETS AND SPACE VOLUMES 


Figure 4-1 provides a functional flowchart for the Earth-refer- 
enced area targets and space volumes. The required inputs are 


a. 

b. 

c. 

d. 

e. 


Target table (either the ground target, ground polygon, or 
ground-fixed space volume table) . 

Target ID. 


Start and end times and time increment (t 
At, respectively) . 

S/C ephemeris and ephemeris ID. 

RNP matrix and associated epoch time. 


start' 


fc end' and 


A brief description of the process is provided herein. The 

heading numbers correspond to the numbered blocks on figure 4-1. 

1. The geodetic parameters defining the area target or space 
volume are obtained from the appropriate target table based 
upon the input target ID. Sections 2.1.2, 2.3.2, and 2.5.2 
describe the specific parameters which are required. 

2. The RNP matrix and its associated epoch time are obtained. 

3. The geodetic coordinates of the area target or space volume 
are transformed to the rotating geocentric coordinate system 
(secs. 2.1.2, 2.3.2, and 2.5.2). 

4. The number of visibility periods between t . t and t gnd 
are determined along with predictions of AOS and LOS times. 
Section 3 describes this process. 

5. A test is performed to determine whether any visibility 
periods were found. If this test is failed, then no AOS/LOS 
t mes are possible and processing is terminated. 

6. A loop is established which will process each visibility 
period found ;tep 4. 
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C«un Ut9*t 



3 


and space volumes. 






7 . 


The current time, 
ion, t. 

zero. The AOS/LOS counter, i, is initialized 


t, is initialized to the AOS time predict- 
'AOS * T ^ e initial value of the precise AOS timer AOS^r 
is set to K 
to one.* 


A loop is established which will terminate when t 
the LOS time prediction, t T __ . 


exceeds 


9. The S/C position vector (at time t in the M50 system) is 
obtained from the ephemeris file based upon the ephemeris ID. 

10. The S/C position vector is transformed to the rotating geo- 
centric system and computations are performed to determine 
whether the S/C lies within the area target or space volume 
(secs. 2.1.2, 2.3.2, 2.3.2, and app. A). 

Based upon the results of these tests, the S/C visibility 
parameter, V, is set 


V > 0 if the S/C lies within the area target or space 
volume 


V < 0 if the S/C lies exterior to the area ..urget or 
space volume 

11. A test is performed on V. 

If V > 0 

11.1 a further test is performed on AOS. to determine 
whether the S/C was also visible during one or more 
of the previous time steps. 

If AOS^ > 0 

11.1.1 then the S/C was visible during one or more 

of the previous time steps. Hence, no transi- 


* Depending upon the target geometry and the S/C groundtrack, 
multiple AOS/LOS times could occur during a single visibility period. 
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tion has occurred. The current time is incre- 
mented by At and the search continues . 

If AOS i £ 0 

11.1.2 then a transition into the area target or 

space volume has occurred. The current AOS 
time, AOS^, is set equal to the current time. 

The current time is incremented by At and 
the search continues. 

If V < 0 

11.2 a further test is performed on AOS^ to determine whether 

the S/C was visible during the previous time step. 

If AOS i > 0 

11.2.1 then the S/C was visible during the previous 
time step and a transition out of the area 
target or space volume has occurred. The 
current LOS time, LOS^ is set equal to the 
time of the previous time step (t - At) . The 
AOS/LOS counter, i, is incremented by one. 

The next AOS time is initialized to zero. The 
current time is incremented by At and the search 
continues. 

11.2.2 the S/C was not visible during the preceding 
time step. Thus no transition has occurred. 

The current time is incremented by At and the 
search continues. 

12. At the completion of the sequential time search, a test is 

made to determine whether the S/C was visible at the end time. 

If the test is true, the last LOS time is set equal to t Tnc . 


4.2 CELESTIAL-FIXED AREA TARGETS AND SPACE VOLUMES 

Figure 4-2 illustrates the logic flow for the celestial- fixed 
area targets and space volumes. This logic is similar to the 
approach presented in section 4.1. The required inputs are 
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a. Target table (either celestial circles, celestial polygons, 
or celestial-fixed space volumes table) . 

b. Target ID. 

c. Start and end times and time increment (t . . , t ., and 

start end 

At, respectively) . 

d. S/C ephemeris and ephemeris ID. 

A brief description of the process is provided herein. The 
heading numbers correspond to the numbered blocks on fiqure 4-2. 

1. The parameters defining the celestial-fixed area target or 
space volume are obtained from the appropriate target table 
based upon the input target ID. Sections 2.2.2, 2.4.2, and 
2.6.2 describe the specific parameters which are required. 

2. The number of visibility periods between t . . and t , 

start end 

are determined along with predictions of AOS and LOS times. 
Section 3 describes this process. 

3. A test is performed to determine whether any visibility periods 
were found. If this test is failed, then no AOS/LOS times 

are possible and processing is terminated. 

4. A loop is established which will process each visibility 
period found in step 2. 

5. The current time, t, is initialized to the AOS time prediction, 
t, n . . The initial value of the precise AOS time, AOS. , is 
set to zero. The AOS/LOS counter, i, is initialized to one. 

6. A loop is established which will terminate when t exceeds 
the LOS time prediction, t 

7. The S/C position vei : u: (at time t in the M50 system) is 
obtained from the v. umeris file based upon ephemeris ID. 

8. Computations are performed to determine whether the °./C lies 
within the area target or space volume (secs. 2.2.2, 2.4.2, 

■^Depending upon the target geometry and the S/C groundtrack, 
multiple AOS/LOS times could occur during a single visibility 
period. 
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2.6.2, and app. A). Based upon the results of these tests, 
the S/C visibility parameter, V, is set 

V > 0 if the S/C lies within the area target or 

space volume 

V < 0 if the S/C lies exterior to the area target 

or space volume 

9. A test is performed on V. 

If V > 0 

9.1 a further test is performed on AOS^ to determine whether 
the S/C was also visible during one or more of the previous 
time steps- 

If AOS i > 0 

9.1.1 then the S/C was visible during one or more of 
the previous time steps. Hence, no transition 
has occurred. The current time is incremented 
by At and the search continues. 

If AOS i < 0 

9.1.2 then a transition into the area target or space 
volume has occurred. The current AOS time, AOS^, 
is set equal to the current time. The current 

time is incremented by At and the search continues. 

If V < 0 

9.2 a further test is performed on AOS^ to determine whether 
the S/C was visible during the previous step. 

If AOS i > 0 

9.2.1 then the S/C was visible during the previous time 
step and a transition out of the area target or 
space volume has occurred. The current LOS time, 

LOS^ , is set equal to the time of the previous 
time step (t - At). The AOS/LOS counter, i, is 
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incremented by one. The next AOS time is initial- 
ized to zero. The current time is incremented by 
At and the search continues. 

9.2.2 the S/C was not visible during the preceding time 
step. Thus, no transition has occurred. The 
current time is incremented by At and the search 
continues . 

10. At the completion of the sequential time search, a test is 
made to determine whether the S/C was visible at the end 
time. If the test is true, the last LOS time is set equal 
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5.0 DETAILED LOGIC FLOW 

Figure 5-1 presents the detailed logic flow for the area targets 
and space volumes processor. This flowchart integrates the 
functional overview presented in section 4 with the equations and 
approach presented in sections 2 and 3 and the appendixes. 

Annotations are provided on the flowchart to describe the comput- 
ations and tests that are performed. In order to simplify 
figure 5-1, several notations have been used to represent corre- 
sponding FORTRAN functions, e.g. , 


<f>(I) 

$(I) 


N (J,I) 


K - 



L = 


1 

1 


represents an array of longitudes indexed by 
the variable I. 

represents a two-dimensional array [i.e., R(3,5)] in 

which the first dimension represents the Cartesian 
components of the vector and the second dimension is 
indexed by the variable I. Hence, operations involving 
vector quantities imply a "DO loop" on the innermost 
dimension. 

represents a three-dimensional array [i.e., N(3,5,3)] 

in which the first dimension represents the Cartesian 
components of the vector and the second and third 
dimensions are indexed by variables J and I, respectively. 
Operations involving these vector quantities also imply 
a DO loop on the innermost dimension, 
implies summation over the indicated range using the 
variable specified by the lower bound (e.g., the variable 
L will be used to perform the summation over the range 
L=1 to L = K - 1) . 
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(5 

Figure 5-1.- Detailed flowchart. 
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Figure 5-1.- Continued 
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Figure 5-1.- Continued. 
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Figure 5-1.- Continued. 
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Figure 5-1.- Continued 
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t A0S (NREV) - t AQSo + (REV-1) T 
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Figure 5-1.- Continued 
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Figure 5-1.- Continued 
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Figure 5-1.- Continued 
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Figure 5-1.- Continued 
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Figure 5-1.- Continued. 
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Figure 5-1.- Continued. 
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Figure 5-1.- Continued 
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APPENDIX A 


SUBDIVIDING CONCAVE POLYGONS 



SUBDIVIDING CONCAVE POLYGONS 


This appendix discusses the procedure for subdividing complex 

concave polygons into two or more simpler convex segments. ^ The 

purpose of this subdivision process is to permit the equations in 

section 2 to be used on a segment-by-segment basis to test for 

containment. The criterion for the S/C (or S/C subsatellite 

point) to be contained in the concave area target or space volume 

is that it must be contained in any one of its segments. For 

convenience, the equations and figures presented in this appendix 

2 

will use a pentagon as an example. However, this approach is 
easily extended to any n-sided polygon. 

Figure A-l presents three examples of concave pentagons. Figure 
A-l(a) illustrates a pentagon having one concave vertex. Figures 
A-l(b) and A-l(c) illustrate pentagons having two concave vertices. 
These pentagons can always be subdivided into triangles by select- 
ing an "appropriate” vertex and connecting nonadjacent vertices 
(fig. A-2). The maximum number of triangles necessary to completely 
subdivide any arbitrarily shaped polygon is 


N 


seg, 


= n - 2 


max 


(A-l) 


whera 


n = number of sides 


A concave polygon is defined to be a polygon that has one or 
more interior vertex angles exceeding 180 degrees. A convex 
polygon is defined to be a polygon that has all of its interior 
angles less than 180 degrees. 

2 

This corresponds to the maximum number of sides specifically 
addressed in the requiremen defined in reference 1. 



(a) One concave vertex. 




(c) Two nonadjacent concave vertices. 
Figure A-l.- Examples of concave pentagons. 


i-2 



Furthermore , the maximum number of interior angles exceeding 180 
degrees can also be determined by noting that the sum of the 
interior angles of the polygon must be equal to the sum of the 
interior angles of all triangles into which it can be subdivided. 
Thus, the sum of the interior vertex angles for any arbitrarily 
shaped polygon is given by 

n 

£ Y i = (n - 2) 180 {hr 2) 

i=l 


where 

= interior vertex angles of the polygon 

Thus, the maximum number of interior vertex angles exceeding 180 
degrees, y*, is given by 

N (Y*) = (n - 3) (A- 3) 

iu»A 

This equation limits the maximum number of concave vertices for a 
pentagon to two. Figures A-l(b) and A-l (c) illustrate two ex- 
amples. In figure A-l (b) , the two concave vertices are adjacent 
to each other. In figure A-l(c), the two concave vertices are non- 
ad jacent. 

The selection of the "appropriate” vertex to begin the subdivision 
process is highly dependent upon the shape ox the polygon and the 
number and relationship of the concave vertices. Also, it is not 
always necessary to subdivide the polygon into triangles. Figure 
A-3 illustrates another method for subdividing the pentagon of 
figure A-lia). In this case, the concave pentagon is subdivided 
into a four-sided convex polygon and one triangle. Furthermore, 
figure A-l by no means exhausts all of the Dotential pentagon 
shapes that could be constructed. 




A- 5 


Since the shape of the area targets and space volumes will remain 
static during a mission, it is recommended that the subdivision 
process be performed manually.* There are two distinct advantages 
to this approach: 

a. It eliminates the coding and execution of complex subdivision 
logic . 

b. It can be performed once for each concave area target and 
space volume and does not have to be repeated each time AOS 
and LOS times are desired. 

The treatment of concave polygons will place additional require- 
ments on the target tables other than those specifically mentioned 
in reference 1. In addition to the number of sides and coordinates 
for each vertex, the target tables must also contain the following 
for each polygon- shaped target 

N - number of segments into which the target is subdivided 
seg 

(3 > N >1 for polygons having five or less sides) 
seg 

VCK - integers defining the counterclockwise ordering of the 
vertices for each segment (i = l,2,...N ge g) 

The use of these additional parameters can best be illustrated by 
example. For figure A-3, this pentagon is subdivided into two 
segments. The first segment is a four-sided polygon defined by 
vertices 1, 2, 3, and 5. The second segment is a triangle defined 
by vertices 3, 4, and 5. The corresponding parameters for this 
pentagon would be 

2 

1235 (or 2351 or 3512 or 5123) 

345 (or 453 or 534) 

*This can easily be performed by plotting the vertex points 
on a Mercator projection. 


N 

seg 

VO, 


vo 2 = 



Similarly, for figures A-2(b) and A-2(c)t 
a. Figure A-2(b) 


b. Figure A-2(c) 


N seg 


3 





V°1 

= 

125 

(or 

251 

or 

512) 

V0 2 


235 

(or 

352 

or 

523) 

V°3 

= 

345 

(or 

453 

or 

534) 

N seg 

= 

3 





v°l 


125 

(or 

251 

or 

512) 

vo 2 


245 

(or 

452 

or 

524) 

V0 3 

= 

234 

(or 

342 

or 

423) 


For consistency, this approach can also be used for convex polygons. 

In this case, N would be one and VO, would be set to the 
seg i 

counterclockwise vertex sequence. 


For computational purposes, the number of sides for each segment, 
n^, can be extracted from the vertex ordering integer, VO^, as 

follows 


= highest values of 


n i 


where TRUNC 


VO 


i 



-1 


>0 


(A-4) 


where 


TRUNC implies integer truncation. 

th 

The vertex numbers, V^ , corresponding to each vertex of the i 

subpolygon can also be extracted from the vertex ordering integer 
as follows 




(A-5a) 


(A- 5b) 
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CONIC INTERSECTIONS 



CONIC INTERSECTIONS 


This appendix presents the equations to compute the intersection 
points between two right circular cones. * The required input 
parameters are 

/V 

- unit vector along the axis of the first cone 

- half-cone angle of the first cone 

A 

- unit vector along the axis of the second cone 
1 2 ~ halftone angle of the second cone 

The quantities to be computed are 

A A 

1^ and I 2 - unit vectors to the intersection points of the 
two cones 


Figure B-l illustrates the two right circular cones and the inter- 
section points. The coordinate system for the cones is arbitrary 

A A 

and the only restriction is that both and must have a 

common origin and must be expressed in the same system. The first 
step is to determine whether the relative orientation of the two 
cones permits any intersections. This can be determined by com- 
puting the angle between the axes of the two cones. 

Y = cos” 1 {V x • V 2 ) (B-l) 


These equations can also be used to compute the intersection 
points between a right circular cone and a pl^ne or the intersection 
points between two planes. 


B-l 




Figure B-l.- Intersection of two cones. 


B-2 



The cones will not intersect if 


Y > y 1 + Y 2 [fig. B-2 (a) ] (B-2) 

or 


Y + smaller {Y^s Y 2 } < greater (Yj/ Y 2 ^ [fig- B-2(b)] (B-3) 

Assuming intersections are possible (i.e., neither equation B-2 
nor B-3 is satisfied) , the next step is to compute the intersection 
vectors. Figure B-3 illustrates the geometry. The intersection 
vectors are symetrically located on either side of the arc connect- 

A A 

ing vectors and V 2 . Furthermore, the intersection vectors 

lie on the perimeter of both cones. Thus, all of the sides of the 

A A A 

spherical triangle connecting vectors V^, V 2 » and are 

known 

A A 

side a (connecting and 1^ ) = Y^ 

A A 

side b (connecting V 2 and 1^ ) = Y 2 

A A 

side c (conni cting V 1 and V 2 ) = y 

The angle from side 'c ' to side 'a ' can be determined using the 
half- angle formula for spherical triangles. The result is 


B = 2 tan 


-1 


K 


sin (S - y 2 ) 


(B-4) 


where 


K 


sin (R - y^) sin(S - y 2 ) sin(S - y) 


sin (S) 


(B-5) 


S = 1/2 (Y x + Y 2 + Y) 


(B-6) 


B-3 




{ 2 



Figure B-2.~ Conditions when cones dc - — 
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B-5 


The dual value of K in equation B-5 produces two values of B 
in equation B-4 . These values of B are equal in magnitude but 
opposite in sign and correspond to the two intersection points.^ - 
The unit vectors to the two intersection points are given by 


where 


A 

1 lx l 2x V lx 


cos B sin 

r i,2 = 

£ ly *2y V ly 


+sin B sin 


£, £ 0 V, 

lz 2z lz_ 


cos y^ 


V 


1 => (V lx’ V ly’ V 


=> U 2x’ l 2y' t 2z ) 

A A 

V, X V 2 
|V X X v 2 \ 


(B-8) 


A 



=> («,, 

lx ly lz 

A /V 

= £ 2 X Vl 


(B-9) 


^For the special case of B = 0, the two intersection vectors 
are coincident. This physically corresponds to the situation when 
the two cones are tangent. 
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